Non-Newtonian thin films with normal stresses: dynamics and spreading. 
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The dynamics of thin films on a horizontal solid substrate is investigated in the case of non- 
Newtonian fluids exhibiting normal stress differences, the rheology of which is strongly non-linear. 
Two coupled equations of evolution for the thickness of the film and the shear rate are proposed 
within the lubrication approximation. This framework is applied to the motion of an advancing 
contact line. The apparent dynamic contact angle is found to depend logarithmically on a lengthscale 
determined solely by the rheological properties of the fluid and the velocity of the contact line. 
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The spreading of a thin fluid layer on a substrate has 
received much attention due to its practical importance. 
However, the motion of a contact line is still a matter of 
debate (see Refs. 0,0,110 for a review). Microscop- 
ically, there is a balance between viscous forces (shear 
viscosity p) and capillary forces (surface tension a). This 
results in the Cox-Voinov law 0, which relates the ap- 
parent (or dynamic) contact angle 9d to the velocity U 
of the contact line 
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x being the distance to the contact line. This equa- 
tion is ill-defined for small x which reflects the diver- 
gence of the viscous stresses at the contac-line 6J. The 
value of the length £ m depends on the regularising mi- 
croscopic physics accounted for in the model — e.g. Van 
der Waals forces 0, slip @ or diffuse interface - 
so that macroscopic measurements can be used to probe 
microscopicproperties. Experiments on the spreading of 
silicon oils Q confirm the model based on Van der Waals 
forces . However it is plausible that the relevant model 
depends on the nature of both the fluid and the substrate. 

In applications, most fluids are complex and exhibit 
non-Newtonian properties. Except for some viscoelastic 
fluids |8|], they have a nonlinear constitutive equation, 
which raises a theoretical challenge 0. Until now, lu- 
brication theories were restricted to fluids with no nor- 
mal stresses, such as yield-stress fluid s llOl or shear- 
thinning fluids [U El El El El El Ellli- Shear- 
thinning was even proposed as the regularising micro- 
scopical mechanism I2I Experimental studies 
are fewer ^1 [51 ; the more recent one E3 a ls° con- 
sidered fluids for which the only non-Newtonian property 
is the existence of normal stresses, for which no theoret- 
ical framework was available. 

In this Letter we consider the spreading of a thin layer 
of fluid having a constant shear viscosity p and exhibiting 
first normal stresses difference [24J a xx — a zz — tp(d z v) 2 , 
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FIG. f: Schematic of the geometry, defining the directions 
x, z, the film thickness h and horizontal velocity v. 



a being the stress tensor (see Fig. 1 for the geometry 
and other notations). In dilute polymeric suspensions, 
the second normal stress difference a yy — a zz is negligi- 
ble and the normal stress coefficient ip can be considered 
as constant |24|. Within the lubrication approximation, 
we propose a set of two coupled equations of evolution 
(Eqs. I11I12|I for the film height and the shear rate (av- 
eraged over the thickness). Then we investigate an ad- 
vancing contact line (at velocity U). In particular we de- 
termine the lengthscale which replaces the microscopic 
length in (JIJ, 
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as a result from the existence of normal stresses, b being 
a numerical constant. Thus we give more grounds to the 
scaling analysis of Ref . E3 ■ 

In view of the lubrication approximation, we introduce 
the aspect ratio e — Z/X of the film, Z being a typical 
thickness and X the horizontal lengthscale. Let p be the 
pressure field and <f> the potential of an applied body force 
such as a microscopic force or gravity (then <j) = pgz, p 
being the volumic mass). The stress balance reads 



dxOxz + d z a zz 



dx4> + dxP 
■ d z (j) + d z p 
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As d x /d z ~ e <C 1, Eq. (0J yields at the lower order in e 
that 7T = p + <p — a zz is a function of x only. Then Eq. 
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(0 becomes d x (a xx — a z 
rheology, 



-d z a zx = d x ii, i.e. , using the 



fid zz v + ip d x [(d z v) 2 ] = d x n. 



(5) 



and 7r is determined using the normal stress balance — p+ 
u zz = jk at the free surface, accounting for its surface 
tension 7 and curvature n, so that 

TT = -yd zz h + (f>(z = h). (6) 

The set ()■")!( it is closed with mass conservation 



d t h + d a 



v(z,t)dz = 0. 
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Thus we obtain a system of PDEs for h and v. In fact, 
Eq. JSJ can be solved for v by a series of the form 



v(x,z,t) = a n (x)z T 



(8) 



where do = to ensure no slip at the substrate, a±(x, t) = 
2s(x, t) is proportional to the mean shear rate across the 
thickness, 02 = d x (—n + 4%/js 2 )/ /z, and each following 
a n can be computed recursively with the ^-derivatives of 
the previous coefficients. Here we propose to truncate the 
expansion at order 2: v = 2sz + a^z 2 . This truncation is 
further discussed in the conclusion, but it can be noted 
beforehand that it is obviously exact in the standard case 
of no normal stress as well as in the case of strong normal 
stress where Eq. shows that the velocity profile is 
linear in z. Then the condition of no shear stress at the 
free surface d z v = yields a second equation relating 02 
to s. As a consequence Eqs. (|SH3> reduce to two coupled 
PDEs for the thickness h and the mean shear rate s, 



2fjts + hd x (ir- 4tps 2 ) = 0, 
d t h + \d x (h 2 s) = 0. 



(9) 
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This set is readily generalised to account for a third di- 
rection y; s = (s x , s y ) is then the vectorial mean shear: 



2^s + /iV(tt-4V>s 2 ) =0, 
d t h+\v- (h 2 s) = 0, 



(11) 
(12) 



where the dynamic pressure n is defined by Eq. and 
V=(d x ,d y ). 

Now we proceed to the study of a moving contact 
line, the only driving forces being the capillary forces: 
7r = —~fd xx h. We consider a contact line advancing at 
constant velocity U towards x — —00, i.e. h{x, t) — 
h(x + Ut) and s(x,t) = s(x + Ut) , so that Eq. fTTTft 
reads U + 2/3 hs = 0. Replacing s in Hm yields 



l + h h2h 
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where C = fxU/j is the capillary number and I = ipU j \x is 
the normal stress characteristic length. The scaling form 
of the solutions to (|13f) is 



x = X£, h(x) =C 1/3 £H(X), 




FIG. 2: Film thickness for an advancing contact line as given 
by the solution to JTSJ. a Thickness H(X). b Slope H'(X) 
(continuous line) and comparison with the asymptotic form 
(91n(6X)) 1/3 with b = 1.69 (dashed line). 



which yields 



1 o ... H 
1 + -H H - 6 — = 0. 
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We look for the solutions of Eq. Ijl5(l which vanish at 
X = and have no macroscopic curvature, i.e. H" — > 
when X — > +00. Such a solution has the expansion 
H{X) = (3/2y/ 3 X 2 / 3 {3+X+aX 2 + 0(X 3 )) near X = 0. 
Shooting on the value of a yields the only solution with 
no macroscopic curvature. It has the classical asymptotic 
form H = 3 2/3 X[\n(bX)}^ 3 , with b = 1.69. This solution 
is depicted in Fig. 2 and allows a matching between Cox- 
Voinov's law (JTJ and the H ~ X 2 ! 3 scaling resulting from 
the balance between capillarity and normal stresses near 
the contact line. 

To summarise, we proposed the set (I11I12|) of coupled 
PDEs for the film thickness and mean shear. It was ob- 
tained using a truncation which is exact in both the lim- 
its of no and strong normal stresses. Within this frame- 
work, we showed that the rheology provides a regularising 
lengthscale £ m (Eq.|2J which is of the order of l^m in ex- 
periments 17]. Obviously the results are valid as long as 



3 



£ m is much larger than any microscopic length such as a 
slip length or the size of a precursor film. Here, the diver- 
gence of the viscous dissipation T> ~ J dX/H is removed 
by the H ~ V 2 / 3 scaling of the film thickness near the 
contact line. The present study could be improved by a 
truncation at higher order, although the robust asymp- 
totics very near to and far from the contact line would 
not be altered; however this would yield a formidable 
numerical task as the order of the PDEs would increase 



with the order of the truncation. Another extension is 
to match the region where normal stresses balance cap- 
illarity to the smaller region where microscopic physics 
become important. 
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